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Abstract 

We analyse the evolution of scalar and gauge fields during a second order phase 
transition using a Langevin equation approach. We show that topological defects formed 
during the phase transition are stable to thermal fluctuations. Our method allows the 
field evolution to be followed throughout the phase transition, for both expanding and 
non-expanding Universes. The results verify the Kibble mechanism for defect formation 
during phase transitions. 
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1 Introduction 



Why there is so httle anti-matter in the Universe, and why the matter coalesced in the way it 
did are two of the major problems facing cosmology. Predictably, both have attracted a great 
deal of attention spawning a panoply of explanations and theories. Some of these theories 
involve objects known as topological defects[|l|, regions of trapped primordial vacuum, an 
example of which is the cosmic string. A string approach is an appealing one since it can be 
used to address both questions: the wake left by strings moving through the Universe can 
produce fluctuations which may lead to the accretion of matter into large scale structures0 0, 
whilst their interaction with particles and the decay of string loops can provide mechanisms 
leading to baryon number violation[Q and the observed matter bias0. Hence, it is important 
to understand how these strings form, both to predict how many we can expect to have been 
created and how likely the above processes are. The formation of topological defects is 
thought to proceed via the Kibble mechanism]^]. 

Modern particle physics and the hot big bang model suggest that as the Universe cooled 
it underwent several phase transitions in which the symmetry of the vacuum was broken into 
a successively smaller, and smaller group. During such a transition, it is possible for fields to 
acquire non-zero vacuum expectation values. How they do this depends on the order of the 
transition. 

If we consider a transition where a U(l) symmetry is broken, then following the transition 
all points in space will have a physically identical, non-zero vacuum expectation value, the 
only variation being in the difference in phase between any two points. 

By causality, we expect the phases to be uncorrelated on distances greater than the 
horizon length, and so there is a finite probability that the phase along a closed path through 
a number of Horizon volumes will wind through some multiple of 27r - an indication that the 
loop contains a string. In practice two points do not need to be separated by a horizon length 
for their phases to be uncorrelated. This should also be true if they are a thermal correlation 
length (defined later) apart; usually a considerably smaller distance than the horizon size. 

This is the Kibble mechanism, and it relies upon the so-called geodesic rule, which is 
that in passing between two domains of different phase, the phase will follow the shortest 
route. In the global case this has been verified both numerically and experimentally!^ P|, 
though, for a local symmetry it has been argued [0] that the presence of gauge fields may 
influence the path the phase takes, and may actually prevent it from following a 'geodesic'. 
More recently, however, work has been done suggesting that, despite this, the geodesic rule 
holds in the local case for a first order phase transition 0. 

As the temperature falls below the critical temperature, for a while it is still possible 
for thermal fluctuations to restore the broken symmetry, and hence erase any topologically 
interesting configuration present at the time. The point at which it is thought that this ceases 
to be possible is referred to as the Ginzburg temperature, T^, and is found by equating 
the free energy with the thermal energy (for such a restoring fluctuation will have a high 



probability while the former is considerably less than the latter). Brandenberger and Davis|ll 



demonstrated that given certain constraints on the parameters, the ratio of fluctuations in the 
scalar field to the background is less than one beneath a temperature just below the Ginzburg 
temperature, regardless of whether gauge fields are, or are not, present. This means that 
topologically non-trivial configurations arising from thermal fluctuations will become stable 
to such fluctuations just under the Ginzburg temperature. This adds weight to the arguments 
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in favour of the Kibble mechanism. 

However, since Brandenberger and Davis considered a hnearized model, only valid for the 
short time immediately following a spinodal decomposition, before non-linear effects start to 
dominate, their analysis only holds for early times. To study the evolution of the fields at 
later times it is necessary to include the non-linearities. A flexible way to do this is to study 
the Langevin equation associated with the classical field equations 1 12]. This is the purpose 
of this paper. 

By studying the Langevin equation for the system, we derive an equation for the prob- 
ability distribution of the fields, P(0j,Af,t) which we use to analyse the evolution of the 
expectation values of the classical fields coupled to a thermal bath. This enables us to study 
not only the stability of configurations to fluctuations at and below the Ginzburg temperature, 
but also the long time evolution of the fields. The flexibility of this method is demonstrated 
by the ease with which it is modified to include the expansion of the Universe. 

Our method seems to be the only one that allows the study of the effect that thermal fluc- 
tuations have on the development of the field, and the stability of defects formed, throughout 
the phase transition. Other methods either concentrate on the start of the phase transition, 
or near its completion. 



2 Global Symmetry 

Although our ultimate aim is to study the case of an expanding universe with broken local 
symmetry, it is beneficial, for several reasons, to start with the more straightforward case 
of a global, non-expanding model. Firstly, since the Kibble mechanism has already been 
verified for this case, we know that any topologically non-trivial configurations present should 
be stable to fluctuations below Tq, and hence we have a benchmark to check our results 
against. Although this provides a useful test of our method, it is by no means a proof of its 
validity. Secondly, since the amount of algebra involved is very dependent on the number 
of fields present, the global, non-expanding case provides the simplest, and hence clearest, 
demonstration of the method used throughout. 
Consider the U{1) toy model 

£ = {D,<P)\Dy) - V{\<P\^) - \f,^F>^^ (1) 

where is a complex scalar field, is the U(l) gauge connection (taken to be zero for now) 
and = — ieA^. We adopt an effective potential of the form 

ni0r) = ^(i0r-^y+ Vi^r, (2) 

where A = (4A + 6e^)/12 and the temperature dependence reflects the fluctuations on a scale 
smaller than some correlation length, defined later. For sufficiently high temperatures this is 
symmetric about a global minimum at zero. However, as the temperature passes through some 
critical temperature, Tq = {X/X)^!], the system undergoes a second order phase transition, 
breaking the U(l) symmetry, with new minima appearing at |0p = rj"^ {1 — T'^/T^). Any two 
points in the new vacuum will now have non-zero vacuum expectation values of equal moduli 
but random phase. 
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Setting = p exp (ia) , our equations of motion for p and a are 

d.d^^p - {d^ad^a)p + = 0, 

(3) 

d^df'a + 29^a9^p/p = 0. 

We assume that a is time independent and varies spatially over some length-scale lnjka- 
Setting R = p we obtain 

+ Fp = , p - i? = 0, (4) 



where 



-\v' + \T^-{k^c.f+'7;{p')\. (5) 




and kc is explained below. Note that we have replaced with A(p^)p/2 c.f the mean 

square approximation to make the resulting equations more accessible. 

For the purpose of this analysis we consider an initial configuration varying spatially 
on a scale of the correlation length. Since we do not want to become embroiled in a dis- 
cussion of effects due to fluctuations on scales shorter than the thermal correlation length, 

= ~ ^g/^c)] ; consider a coarse-grained field where we have integrated 

out all modes associated with such. This leads to the effective potential mentioned ear- 
lier. Hence, if we were to perform a Fourier decomposition, then it would be of the form 
p = J2k<kc Pk {'''k-x) for some k^ ~ We also assume that the mode corresponding 

to kc dominates, (which we later show to be self-consistent) and investigate a configuration 
with a length scale 2-71 /kc- By (j^) we see that the earlier assumption that a = requires 

ka — 2fcc. 

To incorporate thermal fluctuations into our model, we modify such that the equations 
describing the evolution of R and p over some small time interval 6t are 

R{t + 6t) = R{t) -6tFp{t) +6R, (6) 
p{t + 6t) = p{t) +6tR{t) + 6p, (7) 



where 6p and 6R are the thermal fluctuations in p and R respectively. 

Defining pst = pit + 6t), p = pit), and similarly Rst and i?, we can write 

P{psuRsut + 6t) = J d{6p)d{6R)P,{6p)P2{6R) 

xP{pst - 6tRst - 6p, Rst + StFpst - 6R, t) 
d 

— (Pst ~ ^^^st - M 
opst 

d 

x——{Rst + 6tFpst-6R), 



where Pi and P2 are the probability measures for 6p and 6R respectively. 
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Since 5p and 5R are random fluctuations, we may assume that {Sp)~{5R) — 0. Expanding 
the integrands as Taylor series, we find, after considerable algebra, that 



Pst = P-^%(^P^ + ^^M^^PP^ 



d_ 

dp 



d_ 
OR 



(9) 



from which, assuming that 5p and 5R are independent of p and R, we obtain the following 
differential equation for P: 

dP d d 

m = -Op^^'^^SR^'p^^ 

2 5t dp 2 6t OR' ^ ' 



One interpretation of this equation is as follows. If we move the first two terms on the right 
hand side over to the left, then we have a full derivative of P. Liouville's theorem states 
that for a closed system this derivative should be zero. However, our system is coupled to a 
thermal bath and so there is a fiow of probability between the two, as demonstrated by the 
two non-zero noise terms, due to the bath, on the right hand side. 

This equation is clearly a rather forbidding equation to solve analytically. However, we 
can use it to derive equations governing the quadratic moments 

(p^) = / dpdRP{p, R, t)p'^, (pR) = / dpdRP{p, R, t)pR, (i?^) = / dpdRP{p, R, t)R^. 



We choose to investigate these moments because their equations form a closed system with 
the mean field approximation we have taken. 
As a modification we set 

u^{p^)/r^^ , v^{pR)/r^^ , w = {R^) , r^r^t, (11) 

since this normalises u, v and w, and gives the equations for the moments in a form where 
the relative sizes of terms are much more apparent; 

ii = 2v + 5i, (12) 

V = -foU-^Xi^ + w^ (13) 

w = -2foV - Xuv + ^2, (14) 



where ■ denotes ^, 
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are the fluctuations and 



It should be noted that these equations are even nastier than they look at flrst glance, 
since /<, contains a term proportional to T^. However, by assuming that the temperature 
varies at a much slower rate than the flelds (something which we will see is self-consistent 
later on) it is possible to make some progress analytically. 

For the time being, we consider the case when fluctuations are absent, and treat jo as 
constant over a small time period. After a bit of substitution, we integrate the equations to 
obtain 

= -\u^ - ^j^u^ K^u-^ (17) 



where Ai = ^{wi + jo'^i) + Aii^, A2 = 4(t;? — uiWi) and lij, vi and wi are the initial values of 
M, V and w. 

Let the three roots of the polynomial on the right hand side be A*i > /X2 > A*3. Hence 

IIlll2fJ'3^ A*lA*2 + A«2Af3 + A*3A*i = A«i + A*2 + A«3 = (18) 

Taking initial conditions such that, at Tq, Ui = 1 — tc/ta = ^^/[(^(l + \^ / g)], which is the 
minimum of the effective potential at that time, and I'j = = 0, we flnd A2 = 0, implying 
that one of the roots is zero. Furthermore, Ai and /<, are found to be negative, so the two 
non-zero roots must be positive. Since v? is seen to be positive between and ^2 we expect 
u to oscillate between these two. 

We can also make an estimate of the time period of these oscillations, tp, since 

t,= ^A-W, (19) 

where 

//ii /"^i du 



is the flrst complete elliptic integral and k = y [(/xi — 1x2)/ These results closely agree 
with the corresponding numerical calculations. However, they are only valid for very small 
fluctuations, and so we turn to a numerical approach for a more detailed analysis. 

In order to study the effect of ffuctuations on the evolution of the fields, it is necessary 
to make some estimate of the ffuctuation terms. To do this we imagine coupled to some 
other field, ^, in thermal equilibrium, via an extra term in the Lagrangian of the form 

Hi = \gmi^?- (20) 
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By comparing the resulting equations of motion with those already obtained we find that 

Sp = Q , dR = gpip^dt. 



Since ip is in thermal equilibrium, we also have that ip ^ T, with corresponding number 
density = 1.202g^,T^ /n'^ where is the number of internal degrees of freedom (107 for a 
Grand Unified Theory). Taking dt to be a typical interaction time, such that 5t ~ n^^aj^, 
where aj =~ gk~^ is the interaction cross-section, we find the following approximation for 
the fiuctuation terms 



As expected, the size of the fiuctuations decreases with temperature. 

Before we can carry out a calculation, we need to address the problem of choice of pa- 
rameter vahies. Since the model we are considering is a global one. we see from the definition 
of the potential that we must have A = A/3. We now demonstrate why we are justified in 
assuming that the mode corresponding to kc dominates. Figure 1 shows the evolution of the 
quadratic moment u, (corresponding to (p^)) from the Ginzburg time onwards, where we have 
taken A = 0.1, fluctuation coupling, g — A/3 and a range of different wavelengths. The most 
obvious feature is that the mode varying with wavenumber kc dominates those with longer 
wavelengths, consistent with our earlier assumption. However, one may be slightly alarmed 
at the fact that at least two of the curves look like they have no intention of converging to one 
(corresponding to (p^) — > rj'^) as one might expect. The reasons for this are twofold, and both 
somewhat of our own creation. The flrst is that in assuming that p and a vary spatially with 
fixed wavenumber, we alter the value of p for which p = since we have essentially added 
two terms onto the derivative of the potential. This has the effect of raising the equilibrium 
value of p. The second is that we have to make an arbitrary choice of a (taken throughout 
as q; = 1), which effectively scales ka, and so has a similar effect to the first. Conversely it 
could also be used to tune the expected equilibrium value to one by choosing a sufficiently 
small value of a. 

This may raise questions over the validity of this method for studying the evolution of 
fields. However, the evolution of the fields is not qualitatively changed by taking different 
values of kc, k^ or a and so we argue that, as an approximation, our approach is still of 
interest. 

The potential we are using is, unfortunately, only a one loop approximation, and hence 
is not valid above the Ginzburg temperature where higher loops dominate. Our simulation 
therefore must run from the Ginzburg time onwards, and so we can only investigate the 
stabihty of string configurations to thermal fluctuations and not their formation. We also 
only consider the case of a GUT phase transition, since we expect one at a lower temperature, 
such as the Electro-weak phase transition for instance, to be qualitatively the same, but slower 
due to the larger value of tc- We take the coupling A to be between 0.01 and 1, and g < \. 

Figure 2 shows the effect of the two couplings, g and A. The former controls the size 
of fluctuations, and it is seen, in 2a and 2b, that the larger g, the quicker the rise of the 
lower bound, so fluctuations are, bizarrely, actually seen to help stabilise conflgurations, on 








(21) 
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average, by damping oscillations in the field. They also cause the upper bound to rise at 
a faster rate though this is not as pronounced, nor as important to the stability of domain 

structures. 

Fig. 2c reveals that decreasing A decreases the frequency of oscillations, and also the asymp- 
totic value for u. The latter is because kc oc l/^a y/X and we have already noted that 
the value of kc effects the limiting value. Finally, Fig. 2d demonstrates that the effect of 
fluctuations decreases dramatically with A. 

We see then that, since all curves move away from zero, any topologically non-trivial 
configuration is stable from the Ginzburg temperature onwards, though the fields may take 
a long time to reach their equilibrium values. We also note that, in all cases considered, the 
oscillations occur on a much smaller timescale than the evolution towards the equilibrium 
value; consistent with our earlier assumption. 



3 The Effect of Gauge Fields 

That the configurations formed in the above transition are stable against thermal fluctuations 
is nothing new. The Kibble mechanism for the global case is already well accepted, since one 
can argue in favour of the geodesic rule just by demanding that the path followed minimizes 
the energy density. However, the presence of gauge fields may undermine this, since their 
presence in the gradient energy, [Df^(p){D'^(f)y , may make it equally favourable, energetically, 
to follow a longer path. 

Luckily, the method used for studying the global case works equally well in the local one, 
the only drawback being a significant increase in the amount of algebra that has to be done. 
We start by writing the equations of motion in the form 



where q'^ = A'^ — -d^a (Note that this is gauge invariant). Setting Q^^ — q^^, R — p and 
A'' = -dnd'^d'^a for convenience, we find 



where 



R{t + St) = R{t) - StF2{t)p{t) + 5R, 

p{t + St) = p{t) + StR{t) + Sp, 

Q^^(t + 5t) = Q^'lt) - dtG{t)q^{t) + 5Q^' - ^tA^', 

q^{t + 5t) = q^{t) + 5tQ^{t) + 5q\ 



G = {e + 2e^p^)), 



(22) 



and 5R, 5p, 5Q^^, 5q^ are the thermal fluctuations. 
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Proceeding in exactly the same manner as before, we arrive, after some very unpleasant 
algebra, at the equation for P(p, R, g'^, Q^, t). 



2 6t dq^dq'' 2 6t BQ^dQ''' 



It is important to note that the sum over indices in the last two terms is over all four indices, 
not the usual two. Once more we see the violation of Liouville's theorem via the coupling to 
the heat bath. 

From this it is straightforward to obtain the equations governing the quadratic moments. 
Following the global method and defining r, m, f , w as before, plus 

x = {q^)h' , y={q^Q,)h' , z = {Q^)h' (23) 

we arrive at 

ii = 2v + 5i, X = 2y + S3, 

V = — /2M — + e^xu + w, y = — (^) X — 2e^ux + z, (24) 

w = -2f2V - Xuv + 2e^xv + gu (^) ^ + ^2, z = -2 (^) y - Ae'^uy + ^4, 



where 



(k^ A AT2\ 
h = -T-TT + TT-T , (25) 
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and ^2, (^3 and ^4 are the thermal fluctuation terms. Note that the terms involving 
have integrated to zero. Clearly we are not going to get too far with an analytic approach 
this time, so we restrict ourselves to a numerical analysis. 

As before, our first preparation is to calculate the fluctuation terms. Since the coupling 
of (f) to ip is independent of gauge fields, the two new fluctuation terms, ^3 and ^4 must be 
zero. Hence, the only non-zero fluctuation term is ^2, which is unchanged. 

We consider once more an initial domain structure of length scale ^g, and take = 40A/3, 
since we expect the gauge coupling to dominate. Figure 3 shows the effect of varying A and g, 
the results being very similar to those observed in the global symmetry case. Figures 3a and 
3b once again demonstrate the effect of increasing the size of the fluctuations; an increase 
in the rate of growth of the lower bound and a damping of oscillations, whilst Fig. 3c reveals 
how decreasing the size of the self-coupling decreases the frequency of oscillations, the effect 
of fluctuations decreasing in a similar manner (Fig. 3d). 

The most important feature however is that, as in the global case, any non-trivial domain 
structure present at to is seen to be stable against fluctuations at greater times. This rein- 
forces the work by Brandenberger and Davis and, similarly, the arguments in favour of 
the Kibble mechanism. 
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4 Including the Expansion of the Universe 



Until now we have ignored the expansion of the Universe, which we would expect to damp 
the amplitude of any oscillations present. To make the analysis more realistic it is necessary 
to include this expansion. 



4.1 Global Symmetry 



Taking first the case with a global symmetry, such a modification is straightforward, and 
leads to an extra term in the equation of motion for p proportional to the Hubble parameter 



-3H- 



(26) 



dp^ ^ dt 

as expected, a damping term. The only effect this has on the equations for the quadratic 
moments u, v and w, is to add —3Hv to the right hand side of equation (0), and —6Hw to 
that of (p^, though fo is now written as 

2 ^ ^^^^)2 ,_,21 



fo 



A AT2 
2 ^ 2 r/2 



OoV 
a J 



(27) 



where a is the expansion parameter, and oq its value at Tg^. 

Figure 4 shows the results for a small selection of values to illustrate the effects of varying 
the different parameters. Once more we consider an initial domain structure of length scale 
^G- AH four diagrams are seen to display the rapid damping due to the expansion of the 
Universe. 

Figs. 4a and 4b demonstrate the effect of fluctuations. In 4b, where the fluctuations are 
suppressed, the lower bound on u rises much more slowly than in the unsuppressed case. 
Fig. 4a, which actually overshoots its asymptotic value of one before reconverging. Hence, 
it is seen that fluctuations actually make it less likely that a configuration will be erased, 
agreeing with our non-expanding simulations. The upper bound varies very little between 
the two. 

Fig. 4c shows the effects of reducing the self-coupling; a longer period of oscillation, a less 
dramatic initial growth and a much gentler approach toward its asymptotic value. Fig. 4d 
demonstrates how for small values of A the fluctuations have very little effect on the evolution 
of the fields. 

In summary, topologically non-trivial configurations are stable to thermal fluctuations. 
We also note that due to the scale factor now present in the equations of motion, the effect 
of kc and ka is rapidly damped out so that in all expanding cases considered, u converges on 
one, corresponding to {p"^) tending to 1]'^, the long-time minimum of the effective potential. 



4.2 Local Symmetry 

Now including gauge fields once more, the equation of motion for is, predictably, modified 
in a very similar way to that for p when we include expansion, the new version being 

d^d^A" + 2e^p^A'' -2ep^d'a = -^HdoA". (28) 

^ Since a{t) is an unphysical quantity, we can without loss of generality take uq — 1. 
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In addition, the equations for the quadratic moments x, y and z acquire near identical terms 
to those aheady acquired by those for m, v and the only difference being an extra factor 
of two in the former case. 

Once more, we illustrate four different values of the parameters, for a configuration varying 
on a correlation length scale. As in the non-expanding case, we take = 40A/3 throughout. 

Figures 5a and 5b demonstrate the effect of fluctuations. In Fig. 5a, where the fluctuation 
coupling is comparable to the self-coupling, the lower bound to fluctuations is seen to rise 
much quicker (as was the case for a global symmetry) than that in Fig. 5b where the fluctu- 
ations are suppressed. So much so in fact that it overshoots its expected limit, though long 
time studies show that it gradually bends back and converges to one. 

In Fig. 5c we see the effect of decreasing the self-coupling; a less dramatic growth of 
the lower bound, whilst Fig.5d reveals, once more, that the effect of fluctuations decreases 
dramatically with A. Much the same as in the previous three cases. 

Comparing figures 4 and 5, what we notice is that the presence of gauge fields heavily 
damps the initial growth leading to a lower upper bound and consequently smaller oscillations. 



5 Conclusions 

By studying the Langevin equations for the classical fields we have verified that for a U(l) 
model with broken global symmetry (whether expanding or not), string configurations formed 
during a second order phase transition, are stable to thermal fluctuations below the Ginzburg 
temperature. We have also shown this to be true in the case of a system with a local 
symmetry, reinforcing earlier work[|Il|] on the subject, and lending further support to the 
Kibble mechanism for the formation of topological defects. 

The same method has also been used to study how the fields evolve at late times, with the 
scalar field gradually tending to its equilibrium value; a process accelerated by the damping 
produced by an expanding universe. Indeed, our method tracks the evolution of the field 
throughout the phase transition. Other methods are only able to consider early or late times. 
In addition, we have seen that thermal fluctuations actually accelerate the early evolution 
of the fleld, and damp the amplitude of oscillations in the field as it tends to its asymptotic 
value, making it even less likely that a fiuctuation will destroy a configuration. 

This work is still an approximation however, since we have had to assume d = 0, leading 
to an arbitrariness in the asymptotic value of the field in the non-expanding models. However, 
this should not affect the stability of configurations, and in the expanding case this problem 
is smoothed out anyway as the model scales with time. 

Another approximation we have made is in neglecting the dissipation term necessary when 
a source of fluctuations is present |T3| . Since any dissipation term would have a damping 



effect, it would only increase the stability of a non-trivial domain structure, and so including 
it, a further avenue of research, should only strengthen our results. This and the quantum 
field theoretical treatment are in progress ||15| 

Finally, the fiexibility of the Langevin equation approach |]T^ may make this method suit- 



able for a number of other applications, such as the study of the evolution of seed magnetic 
fields following the breaking of a non-Abelian symmetry] IB], and the stability of defects to 
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fluctuations in condensed matter systems such as ^He|T^. Work on these subjects is in 
progress! 
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